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PARACHUTE-DEPLOYMENT-PARAMETER IDENTIFICATION 
BASED ON AN ANALYTICAL SIMULATION 
OF VIKING BLDT AV-4 

By Theodore A, Talay 
Langley Research Center 

SUMMARY 

A six-degree-of- freedom analytical simulation of parachute deployment dynamics 
developed at the Langley Research Center is presented. Voids in the knowledge of vehi- 
cle and decelerator aerodynamics and suspension- system physical properties appear to 
present the major drawbacks to very accurate simulation results. A comparison study 
was made using flight results from the Viking Balloon Launched Decelerator Test (BLDT) 
AV-4. Input parameters unknown or not well defined were iteratively adjusted to force 
an optimum agreement between simulation and flight results of such quantities as vehicle 
velocity, dynamic pressure, Mach number, tension, attitude, and attitude-rate histories. 

On the basis of this study, a set of deployment- parameter input has been defined which 
may be used as a basis for future work in the field of parachute deployment dynamics. 

The results from the study indicate the analytical model is sufficiently sophisticated to 
investigate parachute deployment dynamics with reasonable accuracy. 

INTRODUCTION 

The objective of the Viking '75 mission is to soft-land two science packages on 
Mars. In order to ensure that guidance and control- system capabilities are not exceeded 
during the parachute deceleration phase of the landing mission, limits have been estab- 
lished (ref. 1) for the allowable vehicle attitude and attitude rates. The balloon- launch ed- 
decelerator test (BLDT) program was conducted in the Earth atmosphere for the purpose 
of demonstrating the structural integrity and performance of the Viking decelerator sys- 
tem in a range of flight conditions encompassing those postulated for Mars. These flight 
tests also provided data on the BLDT vehicle attitude and attitude rates. These data are 
useful in determining whether the Viking attitude and attitude -rate limits may be exceeded. 

Although these full-scale tests have been used to demonstrate dynamic performance, 
especially during the parachute deployment phase, analytical deployment models are war- 
ranted since flight conditions at Mars cannot be duplicated precisely on Earth. It is the 
purpose of the present paper to explain one analytical technique, a six- degree- of-freedom 
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model of parachute deployment developed at Langley Research Center4 This technique 
was developed from one designed for stable “inflation parachute studies (ref. 2) and the 
appendix, prepared by the author and W, Douglas Morris of the Langley Research Center, 
defines the modifications made to it to simulate the entire parachute unfurling and infla- 
tion process. 

It will be shown that the model helps define input of vehicle, parachute, and 
suspension- system properties that, for the Viking decelerator system, have been unknown 
up to now or have had high degrees of uncertainty. This is accomplished through the use 
of the model by forcing a good correlation between it and flight-test data from BLDT AV-4. 

Possible future work would include using the entire set of BLDT flight tests to define 
further the system input parameters. Once these inputs were well defined, they would be 
a basis for predicting parachute deployment dynamics and, hence, vehicle attitude and 
attitude rates for Viking missions with Mars flight conditions, 

SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements and cal- 
culations were made in U.S. Customary Units. 
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vehicle aerodynamic axial-force coefficient 
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Cz vehicle aerodynamic normal- force coefficient 

Cz parachute aerodynamic normal-force coefficient 
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— for parachute 
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dx^,dy^,dz. 


damping force in each line, newtons (pounds) 

1 /2 

nominal parachute diameter, , meters (feet) 

diameter of test vehicle, meters (feet) 

column matrix composed of elements dx^, dy^, dz^, meters (feet) 

displacement of suspension- system attachment points with respect to 
initial unstressed position, meters (feet) 



dxo,dyo,dzo 



displacement of the confluence point with respect to initial unstressed 
position, meters (feet) 

force matrix, newtons (pounds) 


jP A^F ^ vehicle aerodynamic forces, newtons (pounds) 

X y z 

F parachute aerodynamic forces, newtons (pounds) 

X y z 

F ^,F ^,F_^ vehicle suspension- system forces, newtons (pounds) 

X y z 

F ^,F ^,F ^ parachute suspension-system forces, newtons (pounds) 

X y z 

Fx^,Fy^,Fz^ forces applied at confluence point, newtons (pounds) 

G ,G^,G 2 gravity components along vehicle axes, meters/s econd^ (feet/second^) 

GxjGyjG^ gravity components along parachute axes, meter s/second^ (feet/second^) 

I^,Iyy,Iz 2 moments of inertia about vehicle axis system, kilogram-meters^ 
(slug-feet2) 

Ixx’Iyy’^zz moments of inertia about parachute axis system, kilogram-meters^ 
(slug-feet2) 

^y’^z’^yz products of inertia about vehicle axis system, kilogram -meters 2 
(slug-feet2) 


i arbitrary suspension-system member 


ij,jj,k 2 unit vector in vehicle axis system 

K spring constant per suspension- system member, newtons/meter (pounds/foot) 

^sec specific secant modulus per suspension- system member, newtons (pounds) 

Kgec specific secant modulus per suspension- system member including effect of 

plastic strain, newtons (pounds) 
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Kg^c D additional specific secant modulus of damping effects per suspension- system 
member, newtons (pounds) 


KrpoT spring constant including effects of plastic deformation and damping per 
suspension- system member, newtons/meter (pounds/foot) 

Krpo'pj matrix of K,pQrji of all suspension-system members 


k: 


TOT 


specific secant modulus including effects of plastic strain and damping per 
suspension-system member, newtons (pounds) 


^o,i 


{^o> 


original unstressed suspension- system member ler^ths, meters (feet) 

matrix of original unstressed suspension-system member ler^ths, meters 
(feet) 


M 


mass of vehicle, kilograms (slugs) 


M 


mass of parachute, kilograms (slugs) 


parachute -system mass components in parachute axis system including that 
due to enclosed and apparent air mass, kilograms (slugs) 

AAA 

vehicle aerodynamic moments, newton-meters (pound-feet) 




A A 
> z 


parachute aerodynamic moments, newton- meters (pound-feet) 


C' 

Mx -“z vehicle suspension- system moments, newton- meters (pound-feet) 

^ ^ ^ ^ 

Mx parachute suspension-system moments, newton-meters (pound-feet) 




vehicle- axis components of vehicle inertial angular velocity, degrees/second 


p,q,r parachute-axis components of parachute inertial angular velocity, 

degrees/ second 


free- stream dynamic pressure, newtons/ mete r^ (pounds/foot^) 
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|Ratejj^ax| absolute value of maximum attitude rate, degrees/second 

Rj inertial position vector of vehicle center of gravity, meters (feet) 

R 2 inertial position vector of parachute center of gravity, meters (feet) 

^Ix I’^ly l^lz 1 scalar components of a vector from one bridle-leg attachment point 

on vehicle to parachute center of gravity in vehicle axis system, 
meters (feet) 

vector from one bridle- leg attachment point to parachute center of gravity in 
vehicle axis system, meters (feet) 

S vehicle reference area, Trd^/4, meters2 (feet^) 

Sq parachute nominal area, meters^ (feet^) 

Sp canopy inlet area, meters^ (feet^) 

T suspension-system tension measured at vehicle, newtons (pounds) 


ijV ^,V 1 X-, Y-, and Z -velocity components of vehicle center of gravity along 

vehicle axes, meters/second (feet/second) 


Vx jjVy ^ X-, Y-, and Z -velocity components of parachute center of gravity 

along vehicle axes, meters/second (feet/second) 

inertial velocity vector of vehicle center of gravity, meters/second 
(feet/second) 

V 2 inertial velocity vector of parachute center of gravity, meters/second 

(feet/second) 

Voo free- stream velocity of center of gravity of vehicle or parachute, 

meter s/second (feet/second) 


X,Y,Z rectangular Cartesian axes 
= Xp^. - - dxQ, meters (feet) 
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AY. = (feet) 

^ ^p,i ■ ^ 0,0 ' meters (feet) 
x,y,z vehicle station coordinates, meters (feet) 

X longitudinal distance of parachute aerodynamic reference center from para- 

chute center of gravity, meters (feet) 


x,y,z X-, y-, and z-distances from vehicle aerodynamic reference to vehicle center 

of gravity, meters (feet) 


X 


Bi-yBi 



inertial velocity of vehicle center of gravity along vehicle axes, 
meters/second (feet/second) 


Xg ,yg ,Zg inertial velocity of parachute center of gravity along parachute axes, 

^ ^ ^ meters/ second (feet/ second) 




inertial wind velocity at parachute in parachute axis system, 
meters/second (feet/second) 




inertial acceleration of vehicle center of gravity along vehicle axes, 
meter s/second^ (feet/second^) 


Xg jYb ’^B inertial acceleration of parachute center of gravity along parachute axes, 

^ ^ ^ meters/second^ (feet/second^) 

Xo Q,y^ q,Zq q original x-, y-, and z-coordinates of confluence point, meters (feet) 

Xp ^,yp ^,Zp ^ coordinates (with respect to parachute axis system) of attachment 

point i including effects of plastic deformation, meters (feet) 


a 


angle of attack (see fig. 2), degrees 


Oij 


total angle of attack (see fig. 2), degrees 


angle of sideslip (see fig. 2), degrees 




elastic deformation of suspension- system member, meters (feet) 
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e elastic strain of suspension- system member, meter s/meter (feet/foot) 

angles between suspension- system member i and parachute X-, Y-, and 
Z-axes, degrees 

,(p. original angles between suspension-system member i in unstressed 
^ ^ ^ position and parachute X-, Y-, and Z-axes, degrees 

/cos ?! - cos ^cos 1//! - cos ^cos < p [ - cos degrees 

r 2 2 2*'l t t » 

i(p’) +('’■') +(<^0 / column matrix composed of the elements pj, t!, and crj 




vector from vehicle center of gravity to vehicle bridle- leg attachment point, 
meters (feet) 


p. p. ^,p- - X-, y-, and z-coordinates of vehicle attachment point of bridle leg 1 

1 xy ^ 1 iz 

in vehicle axis system, meters (feet) 


xi’ yi> zj 


h h 


»^V 

Xi yi zj 


CJ .CO 

"xj’ yf Zi 


inertial vehicle angular velocity components about vehicle axes, 
degrees/ second 

time rate of change of inertial vehicle angular velocity components 
about vehicle axes, degrees/second^ 

inertial parachute angular velocity components about parachute axes, 
degrees/second 

time rate of change of inertial parachute angular velocity components 
about parachute axes, degrees/second^ 


-,6o i,co < angular velocity components of vehicle about vehicle axes, 
x,i y,i z,i / , 

degrees/second 


(JL>- 


angular velocity vector of vehicle, degrees/second 


Subscripts: 


bl 


bridle leg 


10 



fs full scale 

I inertial 

i suspension- system member 

mf mortar firing 

o confluence point 

p parachute 

si stable inflation 

si suspension line 

V vehicle 

Notation: 


first derivative with respect to time 
second derivative with respect to time 
’ includes effects of plastic deformation 

refers to those concerned with parachute 

FLIGHT-TEST DESCRIPTION 

The BLDT system (see fig. 1) consisted of the test vehicle and the Viking decel- 
erator system. The test vehicle was externally similar in configuration to a Viking entry 
vehicle with the exception of rocket motors used to accelerate the test vehicle on powered 
flights. The decelerator system consisted of a 16.15-m (53-ft) nominal-diameter Dq 
disk-gap-band parachute with 48 suspension lines 1.7 Dq in length. Three bridle legs, 
each 2.29 m (7.5 ft) in length, were attached axisymmetrically to the test vehicle. The 
suspension lines and bridle legs met at the confluence point where a swivel was located. 

The parachute, located in the back of the test vehicle, was ejected by a mortar and unfurled 
in a lines-first type of deployment. 
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The BLDT AV-4 test was conducted over the White Sands Missile Range on 
August 13, 1972. The objective of this particular flight was to demonstrate performance 
and structual integrity at deployment conditions in excess of maximum postulated Mars 
dynamic pressure and in excess of Mach 2.0. Actual Mach number was 2.1 and dynamic 
pressure was 498 N/m2 (10.4 Ib/ft^). 

SIMULATION PROGRAM 

The analytical technique used in this simulation is a more general version of one 
originally intended for parachute studies after stable inflation (ref. 2). The new model 
presented here is used to simulate the entire parachute unfurling and inflation processes 
and to predict the resulting influence on the motions of the suspended vehicle. Refer- 
ence 3 describes the model and simulation technique in general and the appendix defines 
the modifications made to the original model from reference 2. 

Briefly, the vehicle and parachute are considered separate point masses constrained 
by a suspension system. Figure 2 indicates the use of massless springs with damping to 
simulate the suspension system. The suspension- line attachment points are constrained 
to move in a plane surface defined by a canopy inlet area. This parameter is input in 
tabular form and simulates the complex flexible nature of the inflating canopy. The vehi- 
cle is considered rigid; thus, the bridle-leg attachment points are fixed with respect to the 
vehicle center of gravity. The juncture of the bridle legs and suspension lines (confluence 
point) seeks a position in space to maintain a static equilibrium of forces. The coordinate 
axis systems are also shown in figure 2. The vehicle and parachute each possess six 
degrees of freedom, that is, changes may occur in the spatial and angular displacements 
in each of the three coordinate directions. In total, therefore, for the entire vehicle- 
parachute system there are twelve degrees of freedom, but it is customary to speak of 
the modeling technique simply as a six-degree-of-freedom technique. 

Input to the simulation consists of masses, moments of inertia, aerodynamic char- 
acteristics and initial state vectors of the vehicle and parachute, suspension- system 
physical properties, and atmospheric data. 

INPUT 

In the ideal case, an analytical model would perfectly describe all the physical 
phenomena that occur and when combined with exact input (i.e., event times, aerodynamic 
coefficients, physical properties, and initial conditions) would precisely simulate the real- 
world event being modeled. In general, however, this is not possible. Complexity, cost, 
and voids in the basic knowledge of decelerator material dynamics and aerodynamic inter- 
actions force the acceptance of a less-than-perfect simulation model and inputs that are 
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suspected of possessing probable errors. When an analytical model is constructed in 
sufficient detail to describe the desired output quantities attention is centered on the 
inputs. An iterative procedure is followed to obtain a good correlation, that is, compare 
desired output quantities with correspondir^ flight data and then revise those suspect 
inputs until a desired degree of comparison occurs. The inputs so derived may not cor- 
respond exactly to those in a perfect analytical model, that is, may not equal exactly the 
real- world-event inputs. Also, there is no guarantee that the accepted inputs are unique; 
there may be other consistent sets of input that yield equally good comparisons. Hope- 
fully, the researcher can justify (based on existing knowledge and experience gained 
through independent comparison simulations) a set as acceptable to the problem at hand. 

Times 

Times for significant events (ref. 4) used in the simulation were: 


Event 

Time from mortar firing, sec 

Initiation of deployment (mortar firing) 

End of unfurling and start of inflation (bagstrip) .... 

Aeroshell separation 

End of simulation 

0.000 

1.233 (ref. 4) 
7.650 (ref. 4) 
50.000 


Mortar firing represents the initiation of parachute deployment and provides the 
reference time for the initial conditions at the beginning of unfurling. Bagstrip represents 
the end of the unfurling process when the simulation switches over to the inflation phase of 
deployment. The time used for bagstrip was derived from flight data. An analytical 
model exists (ref. 5) from which bagstrip times may be determined if conditions at mor- 
tar firing are known. This is of importance when the six-degree-of-freedom model is 
used in a predictive mode in the absence of flight data. At aeroshell separation, the vehi- 
cle physical properties change and the simulation program is reconfigured to reflect this 
at the aeroshell- separation time. The simulation proceeds to 50 sec to ensure that the 
transient deployment motions of the parachute and vehicle are adequately covered. 

Aerodynamic Coefficients 

The aerodynamic coefficients Cy^, Cy^, Cz^, Czp, 

and Cnp for the vehicle and the analogous coefficients for the parachute were unknown 
and assumed equal to zero in the equations for the aerodynamic forces and moments act- 
ing on the system (see the appendix). 

Vehicle.- Other vehicle aerodynamic coefficients for the BLDT AV-4 simulation 
were obtained from wind-tunnel measurements by the manufacturer. The vehicle aerody- 
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namic roll- damping coefficient, not known experimentally, was obtained from iterative 
analysis. A constant value of = -0.125 resulted which gave a close comparison 

between the simulated and flight-test vehicle roll-rate histories. 


Parachute .- The parachute normal-force coefficients and the pitching and yawing 
moment coefficients as functions of the total angle of attack of the parachute were obtained 
from wind-tunnel measurements (ref. 6) at subsonic velocities as no test data were avail- 
able relating supersonic experience. The parachute pitch and yaw damping coefficients 
were not known experimentally. A constant input value of Cmq = Cn^ = -0.20 was 

obtained through iterative analysis, which damped parachute angular attitude rates which, 
in turn, helped match the vehicle attitude and attitude-rate simulations with flight-test 
data. 


The axial-force- coefficient curve for the full-scale parachute was obtained from 
reference 7 but adjusted, through iterative analysis, to higher values in the supersonic 
Mach number range to match the average vehicle -tension- history simulation with flight- 
test values. Additionally, the curve was adjusted at the highest Mach numbers to match 
the first peaks in the tension histories. The results are displayed in figure 3 which pre- 
sents the axial- force coefficients for the full-scale parachute used in this simulation 
along with the curve originally presented in reference 7. From experimental data pre- 
sented in reference 8 for partially inflated parachutes it was determined that the parachute 
axial- force coefficient is, among other factors, a function of the inlet area ratio, Sp^Sp^si. 
A good fit to this functional dependency using the available data was obtained for this sim- 
ulation study using the following relationship: 


Cx 



( 1 ) 


The observed values of atmospheric pressure, density, and speed of sound measured 
just before the AV-4 f Light (ref. 4) were used. Flight-test relative- wind initial conditions 
were input to initiate deployment (at mortar firing). This provided a comparison between 
simulated and flight- test relative -wind vehicle velocity, dynamic pressure, and Mach num- 
ber, However, a no-wind condition was assumed, which meant that relative -wind vehicle 
velocities were assumed equal to the vehicle inertial velocities. The vehicle trajectory 
equations required flight-test vehicle inertial velocities. The simulation substituted the 
relative-wind vehicle velocities instead, and initiated discrepancies in the altitude and 
flight- path angle. The discrepancy in the altitude calculation was updated during the sim- 
ulation to reflect the observed altitude. This,, in turn, ensured that proper values of pres- 
sure, density, and speed of sound were used as they were input as functions of altitude. 

The discrepancy in fUght-path angle was not considered and was believed to have 
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influenced vehicle velocity, dynamic pressure, and Mach number in the latter stages of 
the simulation. A more detailed simulation would include an observed wind input. 

Physical Properties 

Test vehicle .- The test vehicle for the BLDT system is described in detail in ref- 
erence 4. The front of the test vehicle was a blunted 70® half- angle cone with a maximum 
diameter d of 3.51 m (11. 5 ft) and a nominal projected area S of 9.7 m2 (103.9 ft2). 
The vehicle mass, moments and products of inertia, and center-of-gravity location val- 
ues were obtained from measurements by the manufacturer (ref. 4). The physical prop- 
erties of the vehicle at the times an instant after mortar firing and after aeroshell separa- 
tion are presented below: 


Vehicle physical property 

Just after mortar firing 

After aeroshell separation 

Mass, kg (slugs) 

817.3 

(56.0) 

655.3 

(44.9) 

Inertias: 

Ixx, kg-m2 (slug-ft2) . . . 

592 

(437) 

355 

(262) 

lyy, kg-m2 (slug-ft2) . . . 

454 

(335) 

308 

(227) 

Izz, kg-m2 (slug-ft2) . . . 

437 

(322) 

290 

(214) 

Ixy, kg-m2 (slug-ft2) . . , 

0.62 

(0.46) 

0.58 

(0.43) 

Ixz, kg-m2 (slug-ft^) . . , 

2.34 

(1.73) 

5.23 

(3.86) 

lyz, kg-m2 (slug-ft2) . . . 

0.22 

(0.16) 

-0.23 (-0.17) 

Location of center of gravity:^ 
X, cm (in.) 

91.47 (36.01) 

100.71 (39.65) 

y, cm (in.) 

0.00 

(0.0) 

0.00 

(0.0) 

z, cm (in.) 

-3.58 (-1.41) 

-4.47 (-1.76) 


^Zero vehicle station located 8.103 cm (3.19 in.) forward of the virtual vehicle nose 
(ref. 4). 


Parachute.- The parachute for the BLDT system is described in detail in refer- 
ence 4. It was basically a 16.15-m (53.0-ft) nominal diameter disk-gap-band parachute 
with a nominal area Sq of 205.0 m2 (2206,2 ft^). The parachute mass was taken to be 
28.46 kg (1.95 slugs) consisting of the canopy mass and one-third of the suspension- line 
mass. The author of reference 9 suggests adding one-third of the spring mass to the 
mass at the end of the spring to model satisfactorily the dynamic response of a mass 
oscillating at the end of assumed massless springs. In the case of the parachute one- 
third of the suspension- line mass was added to the canopy mass. The center-of-gravity 
locations and moment-of-inertia values for the parachute were computed in the program 
as a function of inputs of canopy inlet area and mass of various components of the canopy 
and as functions of the calculated parachute included and apparent air mass (see ref. 7). 
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Suspension-line properties were obtained in laboratory tensile tests of heat steril- 
ized Type 52, 220-denier, 880-pound minimum tensile strength dacron cord. The tests 

performed at several strain rates ranging from near zero up to 210 percent per sec- 
ond using a single -pull hydraulic ram device. Nonlinear suspension- line elastic and 
damping properties were reported in reference 10. These values were used in the initial 
simulation attempts of the BLDT AV-4 flight test, but it was found that the tensiometer 
time-history output did not compare favorably with flight-test data with respect to both 
the frequency of the oscillations and the envelope of amplitude variations. The 
suspension- line properties were modified by constant factors until the best comparison 
in tension- os dilations frequency and amplitude envelope was obtained. It is believed this 
modification is necessary to simulate accurately the effects of a repeated oscillatory 
dynamic -load history, includir^ time- dependent hysteresis and heating of the material as 
reported in reference 11. These effects would not have been present in the single- pull 
tensile tests conducted in reference 10. From the parametric simulation analysis the 
final values of the suspension- line specific secant modulus per member 1 were 

1.25 times the single -pull values used in reference 7 and the suspension- line damping 
coefficients per member Cgi were 0.75 times the values used in reference 7. These 
properties are presented in figures 4 and 5, respectively, as a function of the suspension- 
line strain. 

It should be noted that 6 members are used to simulate the 48 suspension lines on 
BLDT so that 1 member is equivalent in the simulation to a grouping of 8 suspension lines. 
The values of gj therefore are equivalent to an eight- line grouping of actual sus- 

pension Unes. Also, the values of Cgi are based on an eight-line grouping taken 
together. The values for each BLDT suspension- line specific secant modulus and damp- 
ing coefficients are one- eighth of the values shown in figures 4 and 5, respectively. 

Nonlinear elastic properties for the bridle legs were obtained from the manufacturer 
and the input curve of bridle- leg specific secant modulus per member plotted 

against bridle-leg strain is shown in figure 6. Because bridle- leg damping values were 
unavailable, a viscous damping coefficient per member of C]^]^ = 0 was assumed. This 
is considered a good approximation due to the very large stiffness and the resulting low 
strain and strain rates for the bridle legs. 

Initial Conditions 

At mortar firing for the BLDT AV-4 system the initial conditions for the vehicle and 
decelerator (still in its deployment bag) were equal. The following conditions were 
obtained from references 4 and 12: 
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Vehicle and decelerator parameters 

Value at mortar firing 

Source of data 

Relative velocity, m/sec (ft /sec) . . . 

683 (2 240) 

Ref. 12 

Altitude, m (ft) 

44 862 (147 186) 

Ref. 4 

Flight-path angle, deg 

11.50 

Ref. 12 

Mach number 

2.08 

Ref. 12 

Dynamic pressure, N/m2 (Ib/ft^) ... 

498 (10.4) 

Ref. 12 

Angle of attack, deg 

-4.1 

Ref. 4 

Angle of sideslip, deg 

-3.1 

Ref. 4 

Total angle of attack, deg 

5.2 

Ref. 4 

Pitch rate, deg/sec 

-14 

Ref. 12 

Yaw rate, deg/sec 

4 

Ref. 12 

Roll rate, deg/sec 

-30 

Ref. 12 


Reference 12 was used for some of the inputs since it provided more detailed relative- 
wind histories of the parameters used in the simulation. 

Miscellaneous Conditions 

There are a number of other input quantities which cannot be categorized precisely 
under any of the above headings and are for this reason grouped under the general heading 
of miscellaneous conditions. 

In order to simulate the unfurling process (from mortar firing to bagstrip), the 
analytical model required an input concerning the history of the parachute bag and partially 
inflated canopy drag area C^Sq* This was defined as the vehicle -felt line tension T 
divided by the free- stream dynamic pressure q^ for a massless decelerator system. 

The input history used, obtained by using data derived from reference 4, is presented in 
figure 7. Another required input was the history of suspension- line alinement at the con- 
fluence point relative to the free -stream velocity vector which was grossly approximated 
from data in reference 4. The input values are shown in figure 8. The preceding inputs 
of drag-area and suspension-line-alinement histories during unfurling represented an 
input force and angle-of-force application at the confluence point. From the subsequent 
model analysis the vehicle attitude, attitude rates, and trajectory during the unfurling 
process were output. 

At bagstrip, the analytipal model initiated simulation of parachute motion for the 
inflation phase with the exception of parachute angular rates which were reinitialized. It 
was assumed that the parachute angular rates about the parachute center of gravity were 
equal to the angular motion of the parachute relative to the vehicle, that is, solid-body 
rotation. The values of parachute angular rates about its center of gravity were found, 
from earlier iterative computer simulations, to affect the vehicle angular- rate histories 
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obtained. Therefore, from iterative simulations, the parachute angular rates were 
adjusted to aid in matching the simulated and flight- test vehicle angular -rate histories. 

The values in this simulation were: 

Parachute pitch rate at bagstrip = 15.0 deg/sec 
Parachute yaw rate at bagstrip = 0.0 deg/sec 
Parachute roll rate at bagstrip = -30.0 deg/sec 

The history of canopy-inlet- area ratio throughout the inflation process was a critical 
input parameter. It described the filling sequence and timing as well as breathing motions 
prior to the achievement of a stable canopy. In the simulation of BLDT AV-4, the actual 
inlet- area- ratio history experienced during the flight was used. This inlet- area ratio, 
presented as figure 9, was based on data obtained from onboard movie- camera film, as 
given in references 4 and 13. The inlet-area ratio was corrected to account for the lon- 
gitudinal motions of the canopy inlet resulting from elastic and plastic deformation in the 
suspension lines. A maximum value of plastic deformation equal to 1.8 m (6 ft), based 
on postflight measurements (ref. 4), was used as a best estimate. This was allowed to 
build up as a linear function of time from bagstrip until the loading history stabilized to 
simulate the time -integrated loading effects on the suspension Lines. 

RESULTS AND DISCUSSION 

Due to voids in the knowledge of decelerator and vehicle aerodynamics and physical 
properties, a number of the parameters in this simulation were unknown or had such sig- 
nificant degrees of probable error that they were subjected to iterative analysis. These 
included: 

(1) The vehicle aerodynamic roll-damping coefficient 

(2) The parachute aerodynamic pitch- and yaw-damping coefficients 

(3) The parachute axial-force coefficient 

(4) Suspension-line elasticity 

(5) Suspension- line damping 

(6) Parachute angular rates at bagstrip 

(7) Suspension- line permanent deformation 

The effects of these probable errors are discussed in the following results. 

Figures 10 to 20 present key simulation’ results using the six-degree-of- freedom 
analytical program. Figures 10, 11, and 12 show, respectively, vehicle free- stream 
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velocity, dynamic pressure, and Mach number. Symbols indicate flight-test data and 
uncertainty bands are shown to note tracking and data- reduction inaccuracies. Up to 
about 20 sec after mortar firing the simulated values of free -stream velocity, dynamic 
pressure, and Mach number correspond quite well with flight- test data but thereafter 
there is a notable divergence which is most probably due to the exclusion of winds from 
the analysis. Observed winds (ref. 4) approached nearly 50 m/sec (160 ft/sec) at the 
altitude of parachute deployment and remained nearly so for the time interval of 50 sec 
used in this simulation. The discrepancies in the output become more important as the 
vehicle decelerates and the wind velocities become significant compared to the vehicle 
inertial velocities. 

Figures 13 and 14 show flight- test data and the simulations obtained for the 
suspension-system tension (or vehicle-felt force). Figure 13 represents the suspension- 
system- tension history from mortar firing to 50 sec. After the transient period of major 
canopy breathing (approximately the first 5 sec) the tension histories appear to corre- 
spond quite closely up until about 30 sec after mortar firing after which consistently lower 
tension occurs for the simulation. This is attributed to the smaller simulated dynamic- 
pressure values due to the assumed no-wind condition as discussed earlier. 

Figure 14 gives the suspension-system-tension histories (felt by the vehicle) over 
the first 5 sec after mortar firing where the significant transient behavior occurs. The 
first opening load simulation value of 58 316 N (13 110 lb) agrees closely with the flight- 
test value of 58 952 N (13 253 lb). The tension-history data had a probable error uncer- 
tainty of ±1112 N (±250 lb) for the BLDT AV-4 flight. The very large variations in ten- 
sion and the maximum observed tension value (72 043 N (16 196 lb)) are not as well 
simulated. Several major problems contribute to this result. It is known that a partial 
canopy collapse occurred at about 3,6 sec after mortar firing after a period of relatively 
stable inflation (see fig. 9). This indicates some sort of flow or wake interaction. A 
similar interaction may have caused the very large load variation at about 2.25 sec after 
mortar firing. The parachute axial-force coefficient is a function of a multitude of 
parameters such as Mach number, Reynolds number, body geometry, attitude, wake, and 
body turbulence. For this simulation, however, an average axial- force coefficient as a 
function of Mach number and body geometry was assumed. The detailed effects of 
unsteady flow interactions associated with a breathing canopy and the wake effects of the 
large -diameter test vehicle are not well understood and are not modeled in the simula- 
tion program. Thus, the axial- force coefficient used appears valid for predicting average 
loading behavior but is generally not as good in terms of the large transient-loading 
oscillations. 

As discussed in reference 11, heating and multi cyclic loading and unloading of the 
suspension lines caused hysteresis and plastic deformation. The iterated average values 
of suspension- Line elasticity (fig. 4), damping (fig. 5), and an assumed plastic-deformation 
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history have been used to model these effects. While this may contribute to a realistic 
simulation of the overall nature of the load history in terms of oscillation frequency and 
average loads, a detailed comparison points out the limitations of not having more exact 
experimental data for repeated loadings. For example, up to about 3.7 sec after mortar 
firing, in a region of relatively high suspension- line strain rates, the simulation- loading 
history appears overdamped whereas after 3.7 sec in a region of lower strain rates the 
oscillations appear underdamped. Because of programing difficulties and large experi- 
mental error in the damping data at low strain rates (ref. 10) , the suspension- line damp- 
ing coefficient was assumed to be a function of strain only, and not strain rate. Any future 
simulation in detail would require knowledge of this functional dependency. 

Figures 15, 16, and 17, respectively, show flight-test histories of vehicle pitch, 
yaw, and roll rate and the corresponding simulation histories using the six -degree -of - 
freedom analytical model. The table below shows some major characteristics of the 
pitch and yaw histories to indicate the accuracy of the simulations: 


Vehicle angular -rate characteristic 

Pitch- rate history 

Yaw-rate history 

Flight test 

Simulation 

Flight test 

Simulation 

Time to first node, sec 


9 

15 

17 

Time to second node, sec 


25 

31 

34 

Time to third node, sec 


42 

46 

50 

1 Rate max 1 first node, deg/sec 

111 

105 

82 

82 

|Ratemax| between first and second nodes. 





deg/sec 

46 

58 

34 

38 

|Ratemax| between second and third nodes. 





deg/sec 

28 

38 

26 

27 

|Ft^f6max| between third and fourth nodes. 





deg/sec 

24 

30 

-- 

-- 


Generally the vehicle yaw -rate history is better simulated than the pitch- rate 
history. As indicated in the table, the maximum -rate values calculated up to the first 
nodal point are very close to those observed in flight but the simulated pitch rate damps 
out more slowly than the flight-test data indicates. 

The comparison of flight-test vehicle roll rate and that obtained in the simulation 
shows good results with a close match between the flight- test and calculated nodal points 
(see fig, 17). This indicates that the chosen value for the vehicle roll- rate damping coef- 
ficient (not known experimentally) in this flight simulation is a realistic approximation to 
its true value. It should be noted that the roll- rate histories shown are averaged 
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curves - the high-frequency oscillations have been dropped in the interest of showing 
more clearly the overall long-term nature of the roll- rate history. 

Figures 18, 19, and 20 present the comparison between the flight derived data and 
simulation values of -vehicle a, /3, and Orji, respectively. The flight derived data were 
secured from reference 12. These data were obtained using a trajectory reconstruction 
program to obtain the vehicle Euler angles. These Euler angles were combined with data 
on wind velocity and direction to obtain the values given as flight derived data in fig- 
ures 18, 19, and 20. Generally, the reconstructed values of the Euler angles may be up 
to 5.5° different from cam era- data Euler angles which, tiiemselves, are believed to have 
a mean deviation of ±2° to ±4°. This, coupled with the probable errors in wind data used 
to derive the vehicle a, and leads to some uncertainty in the flight derived data 

presented. 

It is evident that these simulation results do not compare as well as the previously 
presented parameters. Generally, the variations in maximum magnitudes of the angles 
calculated are smaller than for the flight derived values and in the steady state appear to 
assume a regular harmonic oscillatory motion for the angles of attack and sideslip and a 
complex harmonic oscillatory motion for the total angle of attack. Several explanations 
may account for the inaccuracies of the simulation. To begin with, the six-degree-of- 
freedom simulation program appears to give similar general trends and magnitudes in the 
simulation of the vehicle attitude rates but has difficulty in matching flight-test results 
point for point. Integration of vehicle attitude rates gives the angular displacements pre- 
sented. Thus, the errors present in the rate simulation will carry over in the long-term 
integration to angular displacements of the vehicle. 

Inherent to this explanation are the probable errors in the parachute and physical- 
property inputs. Previous computer analysis has given evidence that the parachute motion 
will lead and influence the vehicle motion in a low-frequency mode although the vehicle 
continues to oscillate in attitude at higher natural frequencies commensurate with its 
moments of inertia. The parachute angular rates at bagstrip are important inputs in 
determining subsequent parachute motion and, hence, vehicle motion. The values chosen 
for the inputs are those that appeared to optimize the vehicle attitude- rate comparison in 
overall nature. Their values, however, appear to depend upon conditions of the vehicle at 
mortar firing and would vary from test to test. Further work is indicated to define accu- 
rately parachute motion during the unfurling process. Also, knowledge of parachute aero- 
dynamics is severely limited by a lack of test data - notably aerodynamic force, moment, 
and stability coefficients - performed over a wide range of Mach numbers and possible 
inflation geometries. While available test data were used in the simulation, no values of 
the parachute pitch- and yaw -damping coefficients were available. A value was obtained by 
the iterative process to damp out the parachute motion and, hence, the vehicle motion 
also. Because of the flexible behavior of the suspension system, it is not clear what con- 
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stitutes the parachute’s moments of inertia for pitch- and yaw- computational purposes. 

In this simulation the entire suspension- line moments of inertia were included in the 
computation of the total parachute moments of inertia. These voids in parachute knowl- 
edge lead to undefined uncertainties in the parachute angular-motion simulation and, 
hence, in an indirect sense, in the vehicle angular-motion simulation. 

CONCLUSIONS 

A six-degree -of- freedom analytical model of parachute deployment dynamics was 
developed at the Langley Research Center. A comparison study was conducted using the 
results of the Viking Balloon Launched Decelerator Test (BLDT) AV-4, in an effort to 
define input parameters that, up to now, have been unknown or ill defined. Based on the 
results of this study, the following conclusions are made: 

1. A number of input parameters previously unknown or not well defined, have been 
estimated through the use of iterative analysis. These inputs are considered valid for use- 
as a starting basis for future work involving dynamic analysis of the deployment of a 
supersonic disk-gap-band parachute. In particular, the vehicle aerodynamic roll-damping 
coefficient and suspension-line elasticity are considered to be well defined. The para- 
chute aerodynamic axial- force coefficient, pitch- and yaw- damping coefficients, and the 
suspension- line- damping and plastic- deformation history determined for this simulation 
are considered adequate for use in future simulations. The parachute angular rates, 
suspension-line alinement during unfurling, canopy drag area during unfurling, and the 
canopy-inlet- area ratio during inflation are estimated or known for this particular simula- 
tion. Future work is needed, however, to define these parameters more generally for use 
in predictive analysis, by usii^ test data from the other BLDT flights. 

2. The six-degree-of-freedom analytical model simulates the overall nature of 
supersonic disk-gap- band- parachute deployment dynamics to a reasonable degree of 
accuracy. 

3. Significant voids in the knowledge of decelerator technology, particularly with 
regard to parachute aerodynamic characteristics and suspension- system physical prop- 
erties, appear to be a major obstacle to obtaining very accurate simulations and to the 
use of the analytical model in a predictive mode. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. , July 5, 1974. 
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APPENDIX 


MODIFICATIONS TO STABLE-INFLATION PARACHUTE ANALYTICAL 
TECHNIQUE TO SIMULATE ENTIRE INFLATION PROCESS 

By Theodore A. Talay and W. Douglas Morris 
Langley Research Center 

EQUATIONS OF MOTION 

The equations of motion for a vehicle -decelerator system can be formulated by 
combining the six- degree -of -freedom equations of motion for both the vehicle and para- 
chute subject to the constraint of the suspension system. 

The three translational equations of motion for the center of gravity of the vehicle 
from reference 2 are 
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The three rotational equations of motion of the vehicle about the body axis are 
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The vehicle aerodynamic moments about the center of gravity are 
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The translational and rotational equations of motion of the parachute are different 
from those given in reference 2 since the inflating parachute is nonrigid and has a vari- 
able mass - due to changes in enclosed and apparent air mass (see ref. 7). The three 
translational equations of motion for the center of gravity of the parachute are 
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The parachute aerodynamic forces are 



For the parachute system, the body axes are assumed to lie at the center of gravity of the 
parachute and the axes of symmetry (assuming an axisymmetric parachute) are the prin- 
cipal axes. In this case the products of inertia are zero. The rotational equations of 
motion about the body axes are determined by 



For an aerodynamic reference center lying alor^ the axis of symmetry of the parachute, 
the parachute aerodynamic moments about the center of gravity are given by 
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The parachute-system masses (^,My,M 2 )» mass rates moments of 

inertia (^Ixx>Vy^zz) > moment-of-inertia rates (ixx^yy^zz) calculated using 

a cylindrical-hemisphere model of a variable geometry parachute as given in reference 7. 

It is necessary to determine the suspension- system forces amd moments in order to solve < 
the equations of motion. 


SUSPENSION-SYSTEM FORCES 
Physical Properties 

With the aerodynamic forces and moments known, the problem is finding the forces 
and moments caused by the suspension system. As indicated in reference 2 the axial 
force P in each member i of the suspension system is given by 

Pi = }qei (A9) 

where represents the spring constants and 6^ the elastic deformations of each line. 

If the deformations calculated contain plastic deformations, then a modified value 
of may be found to give the same line forces subject to the same elastic deforma- 
tions or 


Pi = k! 6| (AlO) 

Equating equations (A9) and (AlO) and dividing' the elongations by original suspension- 
system- member lengths ^ to obtain strains, the specific secant modulus that includes 
plastic strain effects is obtained: 
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The damping force in each line is given as 


Di = Ci ^ = qe. 
^0,1 


(A12) 


This may be converted into an additional secant modulus by dividing by c or 


^sec^,i ‘^17 


(A13) 


The final value of specific secant modulus that possesses both an elastic and plastic 
strain effect plus a damping effect is obtained by adding equations (Ail) and (A13) or 
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This may be expressed in terms of a total modified spring constant 
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The values of and an assumed plastic strain history are input. The elastic 

strains Ci are calculated by the model at each time step and are used to call K 

1 seCji 

and from the tabular input where they are stored as functions of elastic strain. 


Attachment Point Distances 

Consider a vector (from the vehicle to the parachute as shown in sketch a) 
which results in the following equation: 

r^ = Rg - Rj - Pj (A16) 
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The vector is fixed in the vehicle (considered a rigid body) from the center of gravity 
to one of the bridle -leg attachment points. Taking the derivative of with respect to 
time gives 

= coj X pj (A18) 


The time rates of change of and R 2 are equal to the inertial velocities of the vehi- 
cle and the parachute centers of gravity: 



The time rate of change of r^ in the vehicle reference frame is 


(A19) 
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Combining equations (A 16) to (A 19) yields 
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If all vector components are expressed in the vehicle axis system there result three 
scalar equations for the time rate of change of the components of the vector from the one 
vehicle bridle -leg attachment point to the parachute center of gravity 
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J 

Equations of the same form are obtained for the time rate of change of the compo- 
nents of the vectors from the other bridle -leg attachment points to the parachute center of 
gravity (expressed in the vehicle axis system). Through transformations., equations (A22) 
are expressible in the parachute axis system. Integration of these equations will give the 
instantaneous positions of the bridle-leg attachment points on the vehicle. 
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The, parachute suspension-line attachment points are also subject to movement with 
time with respect to the parachute axis system as the parachute inflates or breathes. 

The canopy-inlet circle is allowed to expand or contract based upon a history of the 
inlet-area-ratio input. This determines the y- and z-coordinates of the suspension-line 
attachment points with respect to the parachute axis system. Also, based upon the para- 
chute model configuration (ref. 7) durii^ inflation and breathing the inlet circle’s 
x-position with respect to the parachute center of gravity is determined. 

The total deformed lengths of the suspension- system members, expressed in com- 
ponents with respect to the parachute axis system, are given by 


ax! = J 0 - dXo 

= y;,i - yo,0 - dyo > 



- ^ - dz 

0,0 




(A23) 


where the prime (’) values indicate the effects of plastic deformations. Dividing the 
member -component lengths by the total member lengths gives the direction cosines with 
respect to the parachute axis system. 

As shown in reference 2 the force components of each suspension- system member 
with respect to the parachute axis system required to hold the suspension system in 
equilibrium for a given displacement of the vehicle and parachute can be expressed in 
terms of the displacements of the points where the forces are applied (coordinates of the 
suspension- line and bridle-leg attachment points with respect to the original unstressed 
position) and the original lengths and direction cosines of the suspension-system members. 
The matrix equation is given as (from ref. 2) 


0 = [^’][%Ot] [A'f 0 - I (0 0 (t') 2 + (a')2 


(A24) 


where the prime (’) indicates the inclusion of plastic deformations and is the 

matrix of the total modified spring constants determined by equation (A15). 

The confluence point is an interior point in the suspension- system framework and 
it is assumed that no external forces are applied there. Thus, the forces Fx^, Fy^, 
and FZq vanish when the framework is in static equilibrium or 
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(A25) 


The results of combining equations (A15) and (A23) to (A25) are three complex non- 
linear equations (ref. 2) for the three unknowns dx^, dy^, and dz^ - the displacements 
of the confluence point from its original position. The analysis continues as outlined in 
reference 2. A Newton- Raphson iterative solution method is employed to converge on the 
values of dx^, dy^, and dz^. Once the displacements of the confluence point are found, 

the force components at the attachment points for each suspension- system member are 

C C C c 

determined by equations (A24). If these force components ^ 

^ c ~ r ^ 

F , and are multiplied by the moment arms from the vehicle and parachute cen- 
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can be now added to the aerodynamic forces and moments to calculate the six-degree-of- 
freedom motions of both the vehicle and parachute until the next time step. At the new 
time step the position equations (A22) are integrated to obtain new displacements, of the 
bridle-leg attachment points and along with the new suspension- line -attachment- point 
locations obtained from input and parachute -geometry considerations, these are used in 
the determination of the new confluence point position. From this, new suspension-system 
forces and moments are computed and so on. 


The mathematical model as presented has been programed in FORTRAN IV language 
for the CDC 6600 computer. The model is run as an integral part of the LRC-MASS pro- 
gram which is a general multiphase trajectory program. The program starts at mortar 
firing and continues with the unfurling phase until bagstrip when the program changes to 
the inflation phase. The program is later changed again to reinitialize conditions repre- 
senting aeroshell separation. The LRC-MASS program is used to compute the transla- 
tional and rotational motions of the vehicle in three-dimensional space. 
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Figure 1,- Balloon- launched-decelerator test (BLDT) system. 




34 



.2 


.6 


2 . 


Mach number 

Figure 3.- Parachute axial-force coefficient (full scale). 



e, strain 

Figure 4.- Suspension- line secant modulus for eight-line grouping. 
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Figure 6 .- Bridle-leg secant modulus. 
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Figure 7.- Axial-force area during parachute unfurling (input to LRC simulation). 
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Figure 8.- Suspension- Hne attitude during parachute unfur lii^ (input to LRC simulation). 
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Figure 13,- Suspension-system force comparison from mortar firing to 50 sec. 
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Figure 14.- Suspension-system force comparison from mortar firing to 5 sec. 
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Figure 15.- Vehicle pitch-rate comparison. 
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Figure 20.- Vehicle total angle- of- attack comparison. 
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